UNIVERSAL SINGULARITY AT THE CLOSURE OF A GAP 
IN A RANDOM MATRIX THEORY 



E. Brezin 1 and S. Hikami 2 

1 Laboratoire de Physique Theorique, Ecole Normale Superieure, 24 rue Lhomond 75231, Paris 

Cedex 05, France* 
2 Department of Pure and Applied Sciences, University of Tokyo 
Meguro-ku, Komaba, Tokyo 153, Japan 



Abstract 



We consider a Hamiltonian H = Hq + V, in which Hq is a given non- 
random Hermitian matrix, and V is an N x N Hermitian random matrix 
with a Gaussian probability distribution. We had shown before that Dyson's 
universality of the short-range correlations between energy levels holds at 
generic points of the spectrum independently of Hq. We consider here the 
case in which the spectrum of Hq is such that there is a gap in the average 
density of eigenvalues of H which is thus split into two pieces. When the 
spectrum of Hq is tuned so that the gap closes, a new class of universality 
appears for the energy correlations in the vicinity of this singular point. 
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I. INTRODUCTION 



We consider a Hamiltonian which is the sum of a given deterministic part H$ and of 
a random potential V with a Gaussian probability distribution. Although the measure is 
not unitary invariant, one can still obtain the probability distribution for the eigenvalues of 
H through the well-known Itzykson-Zuber integral Jl]. Generalizing a method introduced 
by Kazakov |1 for the density of eigenvalues, we have found an exact representation of the 
n-level correlation functions in terms of the determinant of an n x n matrix whose matrix 
elements are given by a kernel A0v(A, /i) . For generic values of A and \i in the support of 

the average spectrum, we proved earlier that this kernel reduces universally to the sine-kernel 
of Dyson || as if H was not there. Consequently all the correlation functions, including the 
level-spacing distribution, reduce to the Wigner-Dyson form in the short distance regime 
independently of H . However, at singular points of the spectrum the situation is different. 
For instance at the edge of the spectrum of the density of state, the kernel is given in terms 
of Airy functions instead and a new class of universality for the correlations appears • In 
this paper, we investigate what happens when two edge-singularities collapse. The spectrum 
of H is thus tuned to produce a gap in the average density of eigenvalues of H, which closes 
at the origin through a fine tuning of the parameters. The simplest way is to take ±1 for 
the eigenvalues of Hq, with an equal number of positive and negative eigenvalues, but we 
shall prove that the results are independent of H , provided a gap closes. Near the origin, 
a new universal behavior appears, which is not of Airy type. The kernel which governs this 
new singularity will be discussed in detail. A relation to Painleve II differential equations 
and to a two-dimensional A 2 Gamier system is found. Higher multi-critical behavior is also 
investigated. 

We also consider the analogous problem with V made of complex blocks, and in that 
case a similar tuning of H leads to the degeneracy of a Bessel kernel, related to Painleve 
III equations. 



II. DETERMINISTIC PLUS RANDOM HAMILTONIAN 



We consider an N x N Hamiltonian matrix H = H + V, where H is a given, non- 
random Hermitian matrix, and V is a random Gaussian Hermitian matrix. The probability 
distribution P{H) is thus given by 

P(H) = ie-f T - 2 

We are thus dealing with a Gaussian unitary ensemble modified by the external matrix 
source H , which breaks the unitary invariance of the measure. In previous work 0-§J, 
we have discussed the density of state, the two-level correlation function and the n-level 
correlations. For completeness we briefly recall here a few steps, but refer the interested 
reader to our earlier work. The density of state p(A) is 
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p(X) = i < Tr<5(A — H) > 

foe /// 

— e- i7VU f/(t) (2.2) 

oo 27T 

where f/(t) is the average "evolution" operator 

I7(t) =< Tre iNtH > . (2.3) 

We first integrate over the unitary matrix to which diagonalizes H in fl2.1|) ; without loss of 
generality we may assume that Hq is a diagonal matrix with eigenvalues (ai, • • • , a at). This 
is done with the help of the well-known Itzykson-Zuber integral for a unitary matrix u 

f^ M .^M (, 4) 

where A (A) is the Van der Monde determinant constructed with the eigenvalues of A: 

N 

A(A) = H(a i -a J ). (2.5) 

i<j 

We are then led to 

1 N r 

U(t) = „, A ,„, E / dx i • • • dx N e iNtXa A(x 1 , ...,x N ) 

xexp(-^ E^ + ^E ( 2 - 6 ) 
The normalization is fixed by 

17(0) = N (2.7) 

The integration over the Xj's may be done easily, if we note that 

N 



/iV 
dxi • • • dx N A(x!, ■ ■ x N )exp(-— E ^ + ^ E ^ 



A^^-.^^exp^E 6 ?) (2-8) 



Putting 6j = aj + it8 a j, we obtain 



N N I • » „ 

,a a — a™ + zr N ivi 2 



£/(*) = E II ( _ )e-^- +ita " (2.9) 

a=l °7 

The sum over N terms in (|2.9|) may then be replaced by a contour-integral in the complex 
plane, 
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The contour of integration encloses all the eigenvalues a 7 . The Fourier transform with 
respect to t gives the density of state in the presence of an arbitrary external source H ; 
note that this representation is exact for finite N. 



In the large N limit, the one particle Green function G{z) is readily obtained fl2| , |TT 
Indeed one finds that 

G(z) = < — Tr > 

y J N z-Ho + V 

= iV^ z - a 7 - G{z) ^' lV> 

. For simplicity we begin with the case in which half of the eigenvalues a 7 are equal to +a 
and half to —a. When a is larger than one, one finds a gap in the spectrum which closes 
when a reaches one. Then the density of state vanishes at A = from (|2.11|) , p(A) ~ A 1 / 3 . 
The density of state is plotted in Fig.2 of [p| . 
For the n-point correlation function, we have 

Rn(Xi, A 2 , • ■ • , A n ) =< ^Tr5(Xt - H)^Tr5(X 2 - H) ■ ■ ■ ^Tr5(X n - H) > (2.12) 

By using integral representations for the (^-functions, and repeating the previous steps we 
obtain 

R 2 (\i, A 2 ) = K N (X U X 1 )K N (X 2 , A 2 ) - ^(Ai, X 2 )K N (X 2 , A x ) (2.13) 
with the kernel 

^(Ai,A 2 ) = (-l)^ 1 / — /— f[(±L^£ ) _L_ e -fui-ft^mtx 1+ Nux 2 (214) 

J 2tF J 2?!"^ -i U CL^j u tt 

Similarly the n-point functions are given in terms of the determinant of the n x n matrix 
whose elements are K^(Xi,Xj) [p|,[13|. In 0, this kernel K^(Xi, A 2 ) was examined in the 
large N limit, for fixed iV(Ai — A 2 ). In this limit one can evaluate the kernel ( |3.20|) by the 



saddle-point method. The result was found to be, up to a phase factor that we omit here, 

K N (X U A 2 ) = -— sin[7ryp(Ai)] (2.15) 
ny 

where y = N(Xi — A 2 ). Apart from the scale dependence provided by the density of state p, 
the n-point correlation function have thus a universal scaling limit, i.e. independent of the 
deterministic part H of the random Hamiltonian. 

The large-N behaviour of the density of state p(A) = K(X, A) near the edge point is also 
universal. This universality has been investigated for the Airy case and for the Bessel 
case 0. We will show in the next section that similarly at the origin, we find a new class of 
universality for the density of states and for the n-point functions. 
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III. CRITICAL BEHAVIOR NEAR THE ORIGIN 



We first consider the case where the eigenvalues of Hq are ±a, each value being N/2 
times degenerate. The kernel in (|2.14fJ becomes 



iMAi,A 2 ) = (-1)- +1 / ^^(^4)^e-f^f 2 -^^ (3.1) 

J 2tt J 2m u 2 — a 2 u — it 

From this expression, setting Ai = A2, we obtain the density of state p(Ai). The derivative 
of p(A) with respect to A eliminates the factor u — it in the denominator of ( j3.1|) , and leads 
to a factorized expression, 

^P(A) = -0(AMA) (3.2) 

where 

0(A) = r ^e-f* 2 +f H^)-mtx (3 3) 
J-00 2tc 



■0(A) = j> ^e-f « 2 "f K.'-'H*™ (3.4) 

For large iV the two integrals defining the functions <fi and are given by a saddle-point. 
When Ai and A2 are near the origin the saddle-points in the variables t and u move to the 
origin. Therefore for obtaining the large N behavior of <fi near A = we can expand the 
logarithmic term in powers of t. One sees readily that the coefficient of the quadratic term 
in t 2 vanishes for a = 1; in fact three saddle-points are merging at the origin when a reaches 
one. We must then expand in the exponential up to order t 4 and we obtain 

0(A) = *-*«•-*« (3.5) 

J-00 271 

Rescaling t to N~ l l A t' and setting A = N~ 3 ^x we find that 

0( x ) = N 1/4 (f)(N- 3 ^x) (3.6) 
has a large N, finite x, limit given by 
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oc 



Mx) = - I dte~^ tA cos(tx) (3.7) 

7T JO 

This Fourier transform of the exponential of — t 4 is known under the name of a Turrittin or 
Pearcey integral [|14],[15[]. It is immediate to verify that it satisfies the differential equation, 



L'l'i 



x) = x<j){x) (3.8) 



From the integral representation ( |3.7|) we obtain easily the Taylor expansion of this 
function at the origin 
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= V2 - r ( i + f )(-ir2^ 

4* ,£r (2m.)! 



and its asymptotic behavior at large x, 




2 _i _3j ,3^3 4 7T, 



(x) ~ \/ — x-3 e -B xJ cos(— x3 - -) (3.10) 



We return now to the second function ( |3.4| ). Similarly, in the scaling limit, N large, A small, 
iV 3 / 4 A finite, we may expand up to order w 4 , and define 

$(x) = N 1/4 ip(N~ 3 / 4 x). (3.11) 

In the large N, finite x, limit we find 

fa) = I ^A +u *. (3.12) 
Jc 2m 

This integral is no longer a contour integral around the singularities at u — ±1 , but the 
result of a saddle-point approximation. Therefore the integral is over a path consisting of 
four lines of steepest descent in the complex u-plane. Along these straight lines, the variable 
u is changed successively into e "**it and e ± ~ l u. This leads to 



${x) = -Im[- fe"T (e xvw - e~ XUU1 )} (3.13) 

7T JO 

in which uj = eT, The function V>(x) satisfies the differential equation, 

= -xi>(x) (3.14) 
and again we find from fl3.13|) the Taylor expansion 



7 , , 1 ~ (-l) n x 4n+1 (2n)\ 

V> x = = y ± J — ^— J - 3.15 

^to n\(4n + l) K } 



and the large x behavior 



2 ij ,3^3 4 , 2tt x 



^(x) ~ 2 t /— x-ses" 3 cos(^-x3 + y) (3.16) 

In fact one may express the whole kernel K^(\i, X 2 ) of in terms of the two functions 
(f) and ip in the scaling limit. Indeed defining 

\ 1 = N~ 3 ^x, A 2 = N~^y (3.17) 

y) = N^ A K N (N 3 ^Xx, iV 3/4 A 2 ). (3.18) 
in the large N, finite x and y, limit, 
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K(x,y) = — —e * 4 



2tt Jo 2iri 



+ — 7-T - — 7— T) (3-19) 



u — ituo" 1 u — itu u — itcu 3 u — ituj~ 3 
where lu = or more explicitely 



K(x, y) = —- I du j dt 



itx 



2lT 2 i Jo J-oc M 4 + t 4 

x [2m 3 sin a sinh cr + v / 2u 2 t(cos cr sinh a — sin a cosh cr) 
+ 2iut 2 cos cr cosh a — v^2£ 3 (sin cr cosh a + cos cr sinh cr) 



(3.20) 



where a = yu/ y/2. 

In Appendix B, it is shown that the kernel -KT(x, y) may be expressed simply in terms of 
and ip as 

} = HxWiv) - - Ml (3.21) 
x - y 

Therefore the density of state is given by 

p( x ) = -$(x)ip"{x) - 4>"(x)4>'(x) + x4>(x)i>(x)) (3.22) 

and the connected correlation function for two eigenvalues symmetric with respect to the 
origin 

Pc(x,-x) = -\K(x,-x)] 2 (3.23) 



is given through 



K(x, -x) = -4>'(x)il)'(x) (3.24) 
x 



which we know both for x small in an expansion in even powers of x, and for large x, at 
which it behaves as 

2 3 a/3 4 

K(x, -x) ~ - — sm(-^-x^) (3.25) 

OTTX 4 

Using Taylor expansions of 4>{x) and ip(x), we have the expression for the kernel K(x, y) 
by ( p.21|) . In the small x and y, it becomes 

K(x, y) = — -^—\2b + -(-x 2 + 2xy + 2y 2 ) 
2-v^tt 3 / 2 Q 

+ ^(3x 4 - 12x 3 y - 12x 2 y 2 + 8xy 3 - 2y A ) 

+ -^-(-5x 6 + 30x 5 ?/ + 30zV - 40xV + 30xV - 12xy 5 
5040 

-12/) + 0(x s )} (3.26) 
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where a = r(|) = 3.6256 and 6 = = 1.2254. In Appendix C, we evaluate the kernel 
Kj^(x,y) for several finite N cases, and confirmed ( 3.26|) . 

The level spacing probability is an important quantity for the universality. The general 
distribution is given by the determinant of the kernel in our system We consider the 
probability E(s) of no eigenvalues in the interval (— |,f). As same as ordinary random 
matrix theory without an external source, the probability E(s) is given by 



oo / rs/2 f s /2 n 

E ( s ) = E J_ s 2 ■ • ■ j s 2 II dxkdetlKix^x,)}^...^ (3.27) 



f -\Y f s/2 rs/2 
n=0 '" ; "~ s / 2 J- s / 2 k=l 

For the small s, we have a series expansion for E(s) 

rs/2 



E(s) —1 — N f S p{x)dx 

J-s/2 
j\T2 fS /2 fS /2 

+ — / [K(x, x)K{y, y) - K(x, y)K(y, x)]dxdy + ■■■ (3.28) 

2 J-s/2 J-s/2 



Using (|3.12| ), and z = N *A, we have 

E( s ) = 1 - Arl-^r / f (1 + ^Nl z 2 - -N 3 z 4 + ■ ■ Adz + ■■■ 
V ; V2n 3 / 2 J-f 46 8 ; 

= i _ N i-£3L(s + ^iVi(-) 3 - —N 3 (-) 5 ) + ■■■ (3.29) 
v /2 7r 3/2 v ^ 66 V 20 V ; K J 

3 

Then we take a scaling of the energy s as s — > N~^s, 

E(s) = 1 - 42L(s + 4(-) 3 " — (-)') + • ■ ■ (3-30) 
v^vr 3 / 2 66 v 2 ; 20 V 2 ; ; v ; 

Using the expressions of (|3.26|) , we have 

_ — . i \ _ „ c 2 .ab A , a 2 7b 2 . fi 41a6 s . q . .„ „„. 

£(s) = 1 - S + — (— s + ( )s 6 s 8 ) + 0(s 9 ) (3.31) 

K J 2 V 9 V 720 900 ; 201600 ; V ; y J 



where a = T( ±), 6 = 17(f), c = l/(2vr) 3 / 2 and 5 is 

p(x)dx 



rs/2 



s/2 

The quantity s is a function of s. We have to use this variable instead os s since the density 
of state is not a constant in the scaling limit, whereas in the usual Wigner-Dyson case, s = s. 
Therefore, it may be useful to write E(s) as a function of s instead of s. We have in term 
of s, 

E( S ) = 1-~S + ^£!( JL)4 _ (i!^! + Z^)( JL)6 + 0(S 8N 

v ; 18 V 26c ; v 4320 1800 A 26c ; V ; 

= 1 - s + 1.6970s 4 - 16.34555 6 + 0(s 8 ) (3.33) 
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This result may be compared with the Wigner-Dyson case, which reads 

E(s) = l-s + ^-^s e + 0(s 8 ) (3.34) 

The large s behavior may be obtained through a Pade analysis of the series expansion of 
( |3.33j ). In Appendix D, we construct the Pade approximants for the quantity R(s), defined 
by 

R(s) = -^-InE(S) (3.35) 

OS 



and obtain thereby the large s behavior of E(s). 



IV. UNIVERSALITY 

We have analyzed the kernel near the origin when the external source eigenvalues are 
a 7 = ±1. We will show that the results are independent, in the scaling limit, of the distri- 
bution of the external source eigenvalues (provided a gap closes) and are thus universal. We 
consider, for simplicity, the case in which the external source eigenvalues a 7 are distributed 
symmetrically around the origin. If we denote the distribution of a 7 by po( a ) ; the kernel is 
given by 

00 dt r du 1 e _N t 2_N u 2_ Nitzi+NuZ2 



K( Zl ,z 2 ) = - T 



2ixi u — it 



x e f/*wWM$J) (4.!) 
Expanding this expression for small t and u, we have 



x e 

Thus, if we have 



00 dt r du 1 
-oo 2n Jc 2ni u — 

N .,,41 f -J- PO (") 

-(t -u )J dn—r- u<2\ 



K( Zl ,Z 2 ) = - — I — 1 c -^-^ + ^W) f da^l-mtz 1+ Nuz 2 

J-oo 2n J c 2ixi u — it 



da P -^ = 1 (4.3) 
a 2 



we recover the previous kernel ( p.!9|) . The condition of ( |4.3|) is the condition that the density 



of state p(X) starts splitting. Indeed we consider the equation which determines the Green 



function in the large N limit [11,12 



Putting z = 0, noting that parity implies that Re G(0) = 0, when we take the imaginary 



part of (4.4), we obtain 



9 



Therefore, as long as / da^^- < 1 we find that the condition of ( |4.3| ) is equivalent to p(0) = 0, 
a gap is present, but when the condition of ( |4.3| ) is satisfied the gap closes and we recover 
the kernel that we have considered hereabove. Therefore in the scaling regime K(z±, z 2 ) near 
z\ = %2 = is indeed universal, i.e. independent of po( a ) provided the condition (|4.3|) is 
satisfied, or equivalently that a gap closes at the origin in the density of states. 



V. MULTI-CRITICAL BEHAVIOR 

We have obtained the universal behavior of the kernel near the origin by the scaling of 
the energy A = N~ 3 ^x in the large N limit with a fixed x, when the spectrum of H is such 
that the quadratic terms in t in the exponent of the integrand vanish. Now we discuss the 
multi-critical behavior in which the exponent starts from t 2k+2 term in the scaling limit, and 

2fc+l 

the energy A is scaled by A = N 2k + 2 x. 

Suppose that we begin with the Hermitian matrix H coupled to H , which is a diagonal 
complex matrix. The probability distribution of the matrix H is 

Zj 

Then, this probability becomes complex, and we loose the meaning of a probability distri- 
bution since it is not positive definite. We here simply pursue the analytic continuation of 
the diagonal eigenvalues a 7 in the complex plane. By choosing the appropriate a 7 , we can 
obtain an effective kernel, which starts with t 2k+2 term in the exponent of the integrand. 

K ( z _ \ = _ [ ^_ I du 1 -NC(t 2k + 2 +(-l) k u 2k + 2 )-Nitz 1+ Nuz 2 /g 2 ) 

y h ' J 2tt J 2mu-it K ' } 

For k = 2, the eigenvalues of the external source a 7 are assumed to take the values 
a 7 = ±a, ±6, each one beeing ^ times degenerate. For higher multi-criticality, the conditions 
are 

(5-3) 

i< + h ' (5 - 4) 



We then get an effective action for the kernel which starts with a t term ( |5.2| ) in the large 
N limit. The solution of 
form of the kernel is thus 



N limit. The solution of (15.30 and (WM are a = »\ and b = -^e"» % . The effective 



oo 27T / 2Tfi u — it 
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The scaling limit is obtained by rescaling t and u by and A and \i by iV~l. The 

integration of ( |5.5| ) shows a Painleve II equation of A 4 type. 
Similarly to (|3.21|) , we find 

K(x, y) = —I4>""(x)ij(y) - <P"\xW{y) + <P"{xW'{y) ~ <!>'(xW"(y) 
x — y 

+<j>{xW"\y)\ (5.6) 

where <f>(x) and ip(y) are defined by 

^e~- +ltx (5.7) 

-00 Z7T 

Jc 2ni 

The path c of the integration is taken along four lines, starting from the origin to infinity, 
and u is replaced on this lines as u = e ± % l u', and u = e ± ~ l u' . These functions satisfy the 
differential equations, 

d 5 1 

— (j)( x ) = ~-x(f){x) (5.9) 
d 5 1 

_^( x ) = -xifi(x) (5.10) 

In general, if we choose a 7 = (±ai, ±a/j), each cij is N/2k times degenerated, and if 
they satisfy 

k 1 

^ = k 

i=l a t 

£i = 0,(m = 2,3, (5.11) 

i=l a i 

1 2fc+l 

the effective kernel becomes (|5.2|) . The scaling t = N 2k + 2 t , and z = N gives the 

universal behavior in the large N limit. 



VI. DEGENERACY OF BESSEL CASE 

In a previous article |J, we have considered an ensemble of matrices M,of the form 

M=(<i c :) (6.1) 



c 

in which C is an TV x iV complex random matrix, with probability distribution, 
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Noting that 



P(C) = -ex.p(-NTrC^C - NTtAC^C) (6.2) 

Zj 



G(z) = < — Tr — - — > 
y J 2N z-M 

= < -^ T r 2 \^ > (6-3) 



N z 2 - CW 
we obtain a relation between the density of states 



P(z) =< - M ) > (6-4) 

and 

p{z) =< ^Tt5(z - C ] C) > (6.5) 

which reads 

p(z) = \z\p(z 2 ) (6.6) 
The two point correlation function is 

M*U*2) =< ^ Tt6 ^ - ^C)^Tr5(z 2 - C^C) > (6.7) 
Using the Itzykson-Zuber formula, we obtain the kernel in the presence of the external source 

K( Zl ,z 2 ) = r — l ^(— ^)(— ^ TT r 2 ~ U } e- itNz ^ uNz > (6.8) 

When we put a 7 = and take the large N limit, with shifts t,u — > Nt,Nu we obtain the 
Bessel kernel 0. 

If we choose a 1 = — 1 + a, — 1 + b, each one is N/2 times degenerated, and making the 
shifts u — > u — 1 and t — > t + i, the kernel is written by 



dt r du 1 u N (a — it)(b — it) ,n. 



oo 2w J 2tti u — it it (u — a){u — b) 

x e -iNtz 1 +Nuz 2 +N(z 1 -z 2 ) 



dt f du I _ e N M1 .. )+ ^ (1 _»)_» M i. S) _N Ml _i ) 



-oo 2tt J 2-ni u — it 

x e ~iNtz 1 +Nuz 2 +N(z 1 -z 2 ) 

Taking the large N limit with scalings t — > Nt and u —>■ Nu, we have after expanding the 
logarithm to first order, 
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x e -iN 2 tz 1 +N 2 uz 2 +N(z 1 -z 2 ) ]^ 

If a + b is nonvanishing, we have a Bessel kernel, which we have studied in an earlier 
publication. ||. The density of state is obtained as the derivative of K(z,z) ; for example, 
in the case a + b = 2 

dK(z,z) „ /■ dt i h ji f du i ,,„2 

v ' — N 3 — e t' ltx d> — ;e- +,ia; (6.11) 



where we have put z = N~ 2 x 2 . The right hand side of (|6.11 ) becomes N 3 J 2 (2x)/Ax 2 , and 
the density of state in ( |6. 6|) has the universal form, 

p(z) = Nz[J 2 (2zN) + Jl{2zN)} (6.12) 

For the Bessel kernel of (|6.10| ), we have an integral representation. By an appropriate 
rescaling, the Bessel kernel reduces to M 



K(x,y)= <f>(x + zU(y + z)dz (6.13) 
Jo 

where 

*(*) = (6-14) 



Using this representation, and the standard relations between Bessel functions Jq(x) and 
Ji(x), and by considering 



h=l (x + z) JliV¥T " )Jli ^± z) dz (6.15) 



we easily obtain 



h-h=(x- y)K(x,y) 

= ^Ji(Vi) Jo(Vy) - VyJi(Vv)Jo(V^) (e.ie) 

where I2 = I\{x — > y, y — > x). We have thus obtained a closed expression for the kernel from 
this equation; we could have obtained it also through the Christoffel-Darboux identity for 
orthogonal polynomials( Laguerre polynomials in this case). 

When a + b = 0, we find a new scaling limit near the origin. Expanding the logarithmic 
terms in ( |6.9|) , we obtain for the kernel, after rescaling z\ = N~zx 2 , t = N~zt', 

dt' 1 

: exp[- — (a 2 + b 2 ) - itx 2 } (6.17) 



-oo2vr At 2 

This integral should be understood as a contour integral and it should be calculated by taking 
the residue at the origin. This integral / appears in Painleve type III with A2 system, two- 
dimensional Gamier system |17| . The perturbative expansion of this kernel may be deduced 
from this integral representation. 

The higher multi-critical points are obtained by tuning the values of a 7 ; however the 
next one, for instance, requires the conditions; a\ + ai + 0,3 + a 4 = 0, a\ + a| + a| + a\ = 
). Thus, it is necessary to consider the complex values for a 7 , for example, 

d\ — 1, £J2 = —1, CL3 = i, 0-4 = —i. 
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VII. SUMMARY AND DISCUSSIONS 



In this paper, we have considered the critical scaling behavior when the two edges of the 
support of the eigenvalues merge at one point. This has been studied within the Gaussian 
random matrix model with an external source. The correlation functions and the level 
spacing distribution belong then to a new universality class, differnet from the Wigner- 
Dyson case. 

Although our model is limited to Gaussian potentials, several universality classes may 
be obtained by tuning the external source eigenvalues. It may be interesting to extend the 
study to more complicated cases like two dimensional gravity, for which the kernel has a 
similar expression fl8f . 

The unitary matrix model at the transition point has also been shown to be related to 
the Painleve II type by considering non-Gaussian distributions on the unit circle ; multi- 
critical behavior has also been investigated there jnj . Our Gaussian model with an external 
source is an alternative model for tuning multi-critical behavior; it gives concise and exact 
expressions for finite N, which are useful to derive simple closed expressions for the universal 
correlation functions in the scaling limit. 

We have considered the critical case when a gap opens in the density of states. Even 
in this case, the n-level correlation functions are determined by the kernel, completely. 
Therefore, we know that all the correlation functions become universal. After opening the 
gap, there is no more universality. Indeed, if we average over the external source eigenvalues 
in some interval, when the gap opens, we obtain a Poisson like behavior for the level spacing 
probability distribution. 

It has been shown that the Bessel kernel case has interesting physical applications in 
several different problems, including the random flux problem, two-state quantum Hall effect 
flpOfl, and the zero mode of a Dirac operator pifl . It should be interesting to discuss the 



physical applications of the multi-critical behavior discussed in this paper. 
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APPENDIX A: TWO DIMENSIONAL PAINLEVE II EQUATION 



It is known that the differential equation ( |AT| ) is related to a Painleve II differential 



equation with an A2 system, where the notation A2 stands for a two-dimensional Gamier 



system |16| , p!7| . When two saddle points merge, instead of three, the exponent in ( |3.5|) is 
only t 3 , and its Fourier transform is an Airy function; it is then related to Painleve II with 
an A 1 system, and the singurality near the edge of the density of state is indeed described in 
terms of Airy functions [[/]J§. In our case, two edges merge at the origin, three saddle-points 
collapse together, and the Airy case degenerates. 
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Another mathematical interpretation of these functions may be obtained as follows. Let 
us note that the Painleve II equation, is given by the non-linear equation 

^ = 2q 3 + tq + a (Al) 
This equation is equivalent to the two equations, 

|=2p, + a+ i (A3) 

These two equations in turn, may be viewed as Hamilton's equations based on H(p, q) = 
lp2 _ p(g2 _|_ 1^ — -|- i)g. Taking a = — 1/2, p = is a solution; writing then q = 4ln0, 
one finds that <fi satisfies [4>(f)\" = — |0(t) and thus <fi is an Airy function. 

The two-dimensional Hamiltonians for the Painleve II equation (A 2 system) are WJQ 

H\ = (<?2 -Qi- h)pl + 2q 2 piP2 +pj + 2(<?i + t 2 q 2 - 

+ 2( qi q 2 + t x q 2 + t 2 )p 2 + K 9l (A4) 

#2 = q2p\ + 2pip 2 + 2(gig 2 + ^1<?2 + ^2)^1 + 2(<?2 - ?i + ^i)P2 + «<?2 (A5) 
We have, (i, j = 1, 2) 

9tj dpj ' 9tj cty. 
When k = 0, the particular solution is given by 

PiW = (A7) 
Then, w(ti,t 2 ) satisfies the linear partial differential equations, 

d 2 w dw dw 

W = 4tl ^ + 2(2 ^ + 2t " 

d 2 w dw 

4ti7r At 2 w 



2 

If we put t\ = in these equations, we obtain 

9 3 w 



(A6) 



dt\dt 2 dt 2 



8t 2 w (A9) 



In general, if is given by 

w = J^e-^- 2tlu2 - 2t2U du (A10) 
The function in ( ^.8[ ) does satisfy the same differential equation (|A9|) . 
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APPENDIX B: THE KERNEL K(X, Y) 



We have the following equation for K(x,y), which is derived by the contour integral 
representation (|3.1|), 

d ~ * 

—K{x + z,y + z) = -<f){x + z)ifj{y + z) (Bl) 
From this equation, we have 

d * * 

(x - y)—K(x + z,y + z) = -[(x + z) - (y + z)}<p(x + z)^{y + z) 

= -(<£'" (x + z)i){y + z)+ 0(x + z)#'"(y + z)) 
d ~ 

= ~~Q^i ( f ) "( x + z w(y + z ) + <X X + z W {y + z ) - + z )^ {y + z )\ 



(B2) 



Therefore, we get by integration, 



(s — y)K{x + z,y + z) 

= -[<j>"{x + z)$(y + z)+ fa + z)i)"{y + z) - <j>'(x + z)${y + z)\ 

+C(x,y) (B3) 

Setting z = in this equation, it is then easy to prove 

(A + ^) C(l , tf) = (B4) 

since we have 

{-^ + j-)K{x,y) = -4>{x)^{y) (B5) 

Thus, we find from C(x, = C(x — y). Setting further y = in (|B3|), and noting that 
^(0) = V"(0) = and from (pT9| ) 

0) = ft(x)$'(0) (B6) 

we find that the integral constant C(x) is vanishing. 
Thus, we have a simple expression for K(x, y) by 

K(x,y) = (B7) 

x - y 

which is same as ( |3.21| ). 

Since K(x, y) should be finite for y — > x, we have the following constraint between <p(x) 
and ip(x), 
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<j)' (x)ip' (x) — 4>"(x)'if)(x) — 4>(x)ip"(x) = 
This identity is proved by considering the derivative, 

— [(f)' (x)ip' (x) — 4>"{x)^){x) — 4>{x)^)" {x)\ 



(B8) 



-<f/"(x)i/)(x) - <j){x)f"(x) 







(B9) 



Then, by integration, we find the l.h.s of ( P8| ) is a constant, but it vanishes since, by putting 
x = 0, we have 0'(O) = ^(0) = ip"(0) = 0. from @ and (gT5|) . 



APPENDIX C: FINITE N CALCULATION 



We consider some finite N exact calcualtions for the density of state and for the correla- 
tion function near A = 0. 

The derivative of the density of state is factorized as (|3.2| ). We evaluate 0tv(A) of Q3.3|) . 
The expansion at small A is easily done, and after the rescaling x = iV _3//4 A, it gives, for 
instance in the case N = 6, 

(f)(x) = 0(O)(1 - 0.47525a; 2 + 0.09017x 4 + ■ ■ ■) (CI) 



which should be compared with the expansion |3.9| )in the large N limit (f)(x) = 0(0) (1 — 
0.33798x 2 + 0.04166x 4 + - • ■). For we determine the expansion for small A by computing 
the residue in ( |3.4| ) as, 

rtA) = C w jH(-A-i(l)A» + i(l-±)A« 

where A = N^x and CV is a coefficient, which tends to Thus, in the large N limit, 

for a fixed x, we get 

1 2 12 

^(x) = --=x(l--x 4 + -x 8 ---) (C3) 



There are x 3 and x 7 term, but they are order of N~i and neglected in the large N limit. 
This result agrees with fl3.15| ) completely. 

Kn(x, y) is evaluated for small N (N = 2,4,6) by the contour integration (|3.1|) to confirm 
the universal form of (I3.26I). We have 



K 2 (x, y) = ^V^l + (2xy + 2y 2 - x 2 ) + (^ 4 + \xy 3 - 2x 2 y 2 - 2x 3 y + \x 4 ) + • ■ ■] (C4) 
Z7T 6 6 2 
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K ^ y) = ^{\ e ~\ + ^ 2% y + 2 y 2 - * 2 ) 

+ 22& 4 + ^xy 3 - 2x 2 y 2 - 2x 3 y + V) + • • •] (C5) 

66 3 2 

/ \ 1 /7T _q f 87 393 . n 2 2^ 

^e(x,i/) = — ^/-e 3 [- + _(2xt/ + 2y 2 -x 2 ) 

1647, 630 4 4 , 2 2 „ 1 4N . , . 

+ —(j^fy + - W - 2x 3 y + -x 4 ) + • • •] (C6) 

where we have to scale x by N~~ix to compare with (|3.26|) . After the scaling, we find that 
it tends to the universal form in the large N limit. 



APPENDIX D: PADE APPROXIMATION FOR THE LARGE S 



For the Wigner-Dyson case, the probability of no eigenvalue in the interval (— |, |), E(s) 
is Gaussian for large s. Therefore the quantity R(s) defined by 

R(s) = j- g ln E(s) (Dl) 

is proportional to s in the large N limit. 

We apply Pade analysis for R(s) in the Wigner-Dyson case, in which the exact expression 
for R(s) is known. Using the series of E(s) in ( p.34|) , we have 

R(s) = -il + s + s* + (l- ^) S 3 + (1 - ^vr 2 )* 4 + (1 - y + 

7vr 2 7ir\ , . 2tt 2 121tt 4 tt 6 . 7 

+ 1 + s + 1 + s 7 + ... D2 

V 36 675 ; V 9 8100 2205 ; J y J 

The [3,2] Pade approximant for R(s) is given by 

, , 1 + + a 2 s 2 + a 3 s 3 

R{s) = ~ ( i + M + M } (D3) 

where a x = 1.486956, a 2 = 1.90478, a 3 = 0.808162, &i = 0.486956, b 2 = 0.417829. Thus, in 
the large s limit, E(s) is estimated to exp( —0.967097s 2 ). The exact value is known to be 
E(s) ~ exp(-7r 2 s 2 /8), where 7r 2 /8 = 1.233. With the [4,3] Pade, the estimate is further 
improved as a ± = 1.85662, a 2 = 2.56953, a 3 = 1.63848, a 4 = 0.424931, b x = 0.856623, b 2 = 
0.712904, 6 3 = 0.165575, and E(s) ~ exp(-a 4 s 2 /26 3 ), where a 4 /26 3 = 1.2832. This value 
is very close to the exact solution. Using this Pade approximation for R(s), we obtain the 
behavior of E(s) for all region of s by 

E(s) = exp[ [' dxR(x)] (D4) 
Jo 

The level spacing probability p(s) is obtained as the second derivative of E(s). 



18 



p(s) = ^E(s) 

= \^-R(s) + (R(s)) 2 } exp[ f dxR(x)} (D5) 
as Jo 

The value of p(s) by [4,3] Pade analysis for R(s) is shown in Fig.l. 

When we apply this Pade analysis for ( |3.33| ) in terms of s, we have for [3,2] Pade, a± = 
2.1118,02 = 16.0318, a 3 = 9.24379, h = 1.11179, b 2 = 13.920, and E(s) ~ exp(-0.3320s 2 ). 
In the large s limit, s is proportional to s 4 / 3 from ( |3.32|) . The level spacing probability p(s) 
in this case is obtained from the second derivatives of E(s) by s. The result of [3,2] Pade 
for p(s) is shown in Fig.l. 
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